Electron and Spin Transport in the Presence of Complex Absorbing Potential 



Fatih Dogan 1 , Wonkee Kim 1,2 , CM. Blois 1,3 , and F. Marsiglio 1 
1 Department of Physics, University of Alberta, Edmonton, Alberta, Canada, TOG 2J1 
2 Texas Center for Superconductivity, University of Houston, Houston, Texas 77004 
3 Department of Mathematics, University of Toronto, Toronto, Ontario, Canada, M5S 2E6 

(Dated: April 18, 2008) 

We examine the impact of a complex absorbing potential on electron transport, both in the 
continuum and on a lattice. This requires the use of non-Hermitian Hamiltonians; the required 
formalism is briefly outlined. The lattice formulation allows us to study the interesting problem of 
an electron interacting with a stationary spin, and the subsequent time evolution of the electron and 
spin properties as the electron is absorbed after the initial interaction. Remarkably, the properties 
of the localized spin are affected 'at- a- distance' by the interaction of the (now entangled) electron 
with a complex potential. 

PACS numbers: 03.65.Ud,72.25.-b,72.25.Rb,03.67.Mn 



I. INTRODUCTION 

Typically, the dynamics of a closed system are de- 
scribed by a Hamiltonian which is Hermitian and ac- 
counts for all possible interactions and degrees of free- 
dom. However, the task of merely writing down the full 
Hamiltonian is hopeless in many cases and might also be 
superfluous from a practical point of view. Following this 
line of thought, a reduced or effective Hamiltonian can be 
introduced to describe relevant dynamics phenomenolog- 
ically. As discussed in Ref. [l], for example, the effective 
Hamiltonian may or may not be Hermitian, depending 
on the phenomena under consideration. For example, a 
reduction of the full Hamiltonian often used in nuclear 
physics is referred to as the 'optical model'. In order to 
explain experimental results, this model uses phenomeno- 
logical potentials such as a complex absorbing potential. 

When a complex absorbing potential is introduced 
into a system with a real Hamiltonian, the resulting 
effective Hamiltonian is not Hermitian, but symmet- 
ric and complex in its matrix representation. The 
mathematical properties of complex symmetric matri- 
ces are a little more complicated than Hermitian ma- 
trices, as detailed in, for example, Ref. 0. There 
have been many studies, based on various calculational 
schemes, of a complex absorbing potential in the contin- 
uum limit 4ii£i&&£&i&iLi i L 2 . The approach followed in this 
paper requires the diagonalization of a complex symmet- 
ric Hamiltonian defined on a lattice, as will be described 
below. 

The textbook example is the one-dimensional scatter- 
ing problem in which the absorption is modeled with an 
imaginary Dirac-(S potential. The effective Hamiltonian 
is H = p 2 /2m — i/3S(x) with f3 > 0. Then, the scatter- 
ing states are given by (j>k(x) = [e zkx + R e ~ lkx ~\ 0(—x) + 
T e lkx 0(x), where 9 is the Heaviside step function. The 
reflectance and transmittance amplitude R and T are 
determined by matching conditions at i = to be 
R = -p/{2k + (3) and T = 2fc/(2fc + 0). Clearly, 
|i?| 2 + \T\ 2 < 1. In fact, the absorption probability is 
defined as \A\ 2 = 1 — |i?| 2 — \T\ 2 . A deeper understand- 



ing of this 'absorption' can be attained by studying the 
dynamics of a sufficiently broad wave packet; this will be 
illustrated below. 

It is also possible to consider an absorbing potential in 
a (one-dimensional) lattice, by allowing the potential Ui 
at each site i to be complex. This approach follows that 
used for impurity potentials on a lattice as described in 
Ref. Il3l . The reflectance and transmittance amplitude in 
the case of a complex potential U = —i/3 ((3 > 0) at one 
site is, respectively, 
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with the lattice constant set to be unity. These relations 
are the basic building blocks for any more complicated 
complex potential, to be further examined in this paper. 
In particular, spin-flip potentials with the appropriate 
modifications can also be described in this way, to un- 
derstand how spin interactions (along with some deco- 
herence effects) can modify a spin current. The question 
of coherence of spin currents is an important one; in this 
paper we make a preliminary attempt to model decoher- 
ence of spin and charge currents. We find an interesting 
property which highlights the 'entanglement' properties 
of spins, in particular. These issues are the subject of 
this investigation. 

This paper is organized as follows: In Sec. II we present 
the formalism of the dynamics of a wave packet in the 
presence of complex absorbing potentials in both the con- 
tinuum limit and on a lattice. Sec. III. illustrates how 
a complex potential influences the time-evolution of a 
wave packet as a simple example. In Sec. IV we study 
the spin-dependent scattering problem in the presence of 
a complex absorbing potential. The time evolution of 
an expectation value of an operator is derived when the 
Hamiltonian is non-Hermitian. We then demonstrate a 
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highly non-trivial consequence of the complex potential 
on spin transfer. This new result shows that quantum 
entanglement plays an important role for spin transfer 
in the presence of a complex potential. In Sec. V we 
conclude with a summary. 



II. FORMALISM 

We begin by considering a single particle in d- 
dimensional space under the influence of a complex po- 
tential. As usual, the state space is L 2 (R d ) equipped 
with the usual inner product 



(/Is) == / f*g 



(3) 



The Hamiltonian, H := p 2 /2m + V where V : R d -> C 
achieves non-real values (on some set of non-zero mea- 
sure), is easily seen to be non-symmetric (and therefore 
non-Hcrmitian) . Computational complications due to 
a non-Hcrmitian Hamiltonian are many; some nuances 
have been catalogued in Ref. [T3 and Ref. d, for exam- 
ple. We follow Moiseyev^ and introduce the c-product, 
(■!■): L 2 (R d ) x L 2 (W l ) -> C, defined by 



(/Iff) := / fg 



(4) 



(Note that by the Cauchy-Schwarz inequality, |(/|g)| = 

K/*|S)I < II/1IIMI = ll/IIIMI < °° for f,g e L 2 (R d ).) 
We see that the Hamiltonian exhibits a symmetry with 
respect to the c-product; that is, for all functions / and 
g in the domain of H (the Sobolev space if 2 (R )), 



(f\Hg) = (Hf\g). 



(5) 



This can be easily shown for / S H 2 (M. d ) and g S 
C£°(R d ), and a density argument completes the proof. 
From this symmetry, it follows that eigenfunctions of 
H corresponding to distinct eigenvalues are 'orthogonal' 
with respect to the c-product (rather than the inner prod- 
uct [3] as usual). That is, if ip, (f> € L 2 (R d ), A, £ € C with 
A 7^ £, and Hip = Xip, H<fi — then 



(V#) = o. 



(6) 



For convenience, we now restrict our domain to a d- 
dimensional cube in R d and impose periodic boundary 
conditions on the edges. In doing so, we force the spec- 
trum of H to be discrete. We look for solutions to the 
initial value problem, 



(7) 



We suppose that, as in the case where H is Hermitian, the 
solution is obtained by acting on the initial wave function 
(f> with the time-evolution operator U(t) := er lHt . Then 



if the eigenfunctions of H are denoted {4>j}°^ 1 with cor- 
responding complex eigenvalues Ej , and the initial wave 
function is ip = ^ a 
Eq. © is 



3 ' 

it follows that the solution to 



-iEjt 



3 ■ 



(8) 



See Appendix A for a sample calculation of wave-packet 
evolution under the influence of a complex potential. 

In a lattice with N sites, we use the tight-binding 
Hamiltonian: 



E — -t ^ C ia G ic 
<i,j>a 



(9) 



where Cf^ creates an electron with a spin a at a site i, 
to is a hopping amplitude between the nearest neighbor 
sites and is set to be unity, Uj is a complex potential at 
site j, and X represents a set of the potential sites. For 
complex absorbing potentials, Uj — —i/3j, with j3j > 0. 
Varying Ui and X, one can build various effective poten- 
tials. For example, in order to study the dynamics of a 
wave packet and compare the reflectance [Eq. (JTJ)] and 
the transmittance [Eq. (|2"])] given in the Introduction, we 
choose a single-site potential with Uj = —if3So.j- 

In the lattice, the state space is finite-dimensional; the 
Hamiltonian can be represented by a complex symmetric 
matrix. Recall that an N x N matrix with entries in a 
field (C) is diagonalisable if and only if it has N linearly 
independent eigenvectors. 

An initial electron wave packet with spin up can be 
written as 



/V 



|*(0)H]>>^l ) 



(10) 



where Ci acts on the vacuum |0) to create an electron 
with spin up, and |<^i| 2 is the initial probability to find 
such an electron at the site i. Thus in the matrix repre- 
sentation 



l*(0)> = 



/ Vi (0) \ 
Va(0) 



V Mo) / 



(ii) 



For the purposes of this paper, 1^(0)) is a Gaussian wave 
packet which is far away from (and broad when com- 
pared to) the region of non-zero potential. In order to 
find the time evolution of this wave packet, we solve for 
all eigenvectors |n) and corresponding eigenvalues E n of 
the Hamiltonian matrix (and hence verify its diagonalis- 
ability). Now the time-evolution of the wave packet on 
the lattice is given by 



Wt)) 



N 



(n|*(0))e" li5 " t |n 



(12) 
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where (n| is the matrix transpose of |n). This can be seen 
by noting that the eigenvectors \n) form an orthogonal 
set with respect to the finite-dimensional analogue of the 
c-product: 

(n\m) = S mn (13) 

for all m,n £ {0, N}. 

The reflection and transmission probabilities are de- 
fined as \TZ\ 2 := lim^ Ei<mini l<ft(*)| 2 and l T l 2 : = 
limt^oo J2i>maxi lw(*)| 2 respectively. Here, ^(t) is the 



:th 



component of \^(t)), and X is the region of non-zero 
potential. 

So far, we have acquired the machinery to describe the 
dynamics of a given initial wave function in the presence 
of a complex potential by solving the Schrodinger equa- 
tion. In the next sections we will utilize it to examine 
various physical systems. 
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III. COMPLEX POTENTIALS AS ABSORBERS 

For illustrative purposes, a good example is the com- 
plex Dirac-5 potential, V(x) = XS(x), where A <E C with 
Im(A) < for an absorbing potential. (The complex 
square-well potential is also described in detail in Ap- 
pendix A.) Suppose that the one-dimensional space spans 
from x = —L to x = L. Then, eigenfunctions of the 
Hamiltonian, H = p 2 /2m + V(x), can be classified ac- 
cording to their parity. These eigenfunctions must have 
the form 

J Ccos(fca;) + Dsin(fca;) if — L < x < 
^ ~ 1 A cos(fcr) + B sm{kx) if < x < L] , ^ ' 



where, if E is the corresponding eigenvalue, k = \j2mE. 
For even eigenfunctions, C — A and D = —B. We 
impose open boundary conditions at x = ±L (namely, 
ip(—L) = ip(L) = 0). Using the usual matching condi- 
tions for the eigenstates at x — 0, we find that 



, s 2k 
taii(kL) = — — . 

A 



(15) 



For the odd channel, C — —A, D = B, and the eigen- 
states should vanish at x = 0, which requires that A = 0. 
As expected, the odd-channel eigenstates are indepen- 
dent of A. These states are given by ip(x) — sin(p n x) 
with p n — nir/L. 

The initial wave packet, 



*o(a:) = (2?ra) 



-l/4 g ifc (x+a;o) e 



(16) 



can then be expanded in terms of the eigenfunctions of 
H by exploiting their c-product orthogonality: = 
J2k c ki>k, where 



Cfe 



W*l*o). 

The time evolution is then given by Eq. 



(17) 



FIG. 1: (Color online) The time evolution of a wave packet 
in the presence of a complex absorbing potential. 



For the lattice problem, as described in Sec. II, one 
can consider a complex potential with U — —i(3 at 
a single site, say /. Then the Hamiltonian becomes 
H = -toJ2< ilj>r ,C+C ja - iPZ.CtC Ia . Since the 
Hamiltonian is a complex matrix, the numerical diagonal- 
ization can be done by the driver ZGEEVX contained in 
the LAPACK packaged In Fig. 1, we show the time evo- 
lution of a wave packet, which initially resides at xq = 50 
with the average momentum ko = ir/2. The single im- 
purity is at / = 100 with (3 — 1. Fig. 2 compares the re- 
flectance and transmittance of the wave packet with the 
results from the scattering plane-wave approach for var- 
ious (average) momenta ka. The scattering plane- wave 
approach is, in fact, verified by the dynamical calcula- 
tions. As long as the wave packet is broad enough, the 
scattering plane- wave approach for a complex square well 
can be verified in a similar manner. 



IV. EFFECTS OF A COMPLEX ABSORBING 
POTENTIAL ON SPIN TRANSFER 

We now examine the interaction of an electron with a 
spin-flip potential followed by a complex absorbing po- 
tential. In the continuum, a simple representative Hamil- 
tonian is H = p 2 /2m — 2Jqij ■ SS(x) — if35{x — a), where 
Jo is the coupling constant of a spin-flip interaction be- 
tween the electron spin a and a local spin S at x = 0. The 
distance between the two interactions, a, is taken to be 
much larger than the width of an incoming wave packet. 
Note the conflicting requirements; a broad wave packet 
is required to have a reasonably well-defined momentum, 
but a narrow wave packet (on the scale of a) is required 
to observe causal effects properly. The wave packet inter- 
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FIG. 2: (Color online) The reflectance and transmittance 
probability as a function of the momentum. The scattering 
approach is verified by the full dynamics. 



acts with the local spin first and is partially reflected and 
partially transmitted. (Note that the amplitudes of the 
transmitted and reflected waves depend upon the initial 
spin configuration.) The transmitted wave later interacts 
with the complex absorbing potential after completely 



emerging from the spin-flip potential. (In other words, 
the wave packet does not simultaneously feel both inter- 
actions at any given moment.) 

It is numerically feasible to find some eigenvalues 
for this Hamiltonian; however the comprehensive search 
for all complex eigenvalues within some upper modulus 
bound (given by appropriate truncation limits) is seen 
to be a formidable task. The analogous problem on the 
lattice is far more tractable, owing to the fact that the 
Hamiltonian is a complex symmetric matrix, which is di- 
agonalizable by standard numerical routines. The lattice 
Hamiltonian is 

H = —to CfoCja — 2J ai ■ S; + 2^ UiC^G'io- , 

<i.j>(7 2(zX,0- 

(18) 

where the local spin is at the site I and the complex 
absorbing potential sites belong to the set I far away 
from the site I. Since we have only one local spin on the 
lattice, we will drop the site label I from now on. 

We will investigate the case where initially the local 
spin is pointing down (S z = —S) while the electron spin is 
pointing up (er 2 = 1/2). In particular, we wish to monitor 
the time evolution of (S z ) as the electron interacts with 
the local spin and is subsequently partially absorbed by 
the complex potential. 

The expectation value of a time-independent operator 
A is (A) = (y(t)\A\V(t))/{y(t)\V(t)}. Differentiating 
with respect to time, 



d (*\A\*) _ (*|*)4(*|>t|*)-(*|>t|*)4<*|* 



dt (#|#) 

[H+A 



AH] |*) 



-i{^\A\^){^\ [H+ 



mm 



i(H + A - AH) 



(\Br|¥)2 

i(A)(H + - H) 



(19) 



Note that (A) is not guaranteed to be constant in time, 
even if [H, A] =0. If H is Hcrmitian, this equation re- 
duces to the usual expression: d(A)/dt = i([H,A\). For 
example, the z-component J z of the total spin J = a + S 
commutes with the Hamiltonian: [H, J z ] = 0. However, 
it is clear from Eq. (fl9|) that in general ^{J z ) ^ 0. Con- 
sequently, in general, we cannot utilize (J z ) as a con- 
served quantity. We will further discuss this issue below. 

The spin space of this system has a basis of 2(2S f + 1) 
spin states: \a z ,S z ) = | + , S) 7 ■ ■ ■ , | — , — S). We focus on 
the simplest case, S ~ 1/2, in which there are four spin 
basis states: |+,|), | + ,|), |-,t), The Hamilto- 

nian matrix is of dimension (4iV x 47V), where N is the 
number of lattice sites. The initial electron wave packet 
is also expressed in terms of these basis states; 



l*(o)> = X>u(o)|+,T> + ^(o)|+,|) 

3 

+¥>3j(0)|- T> +^-(o)|-i) (20) 

For most of this investigation, only if2,j(0) is non-zero, 
with 

y, 2J (0) = e ' fc °(j-J°) e ~°4i o)2 i (21) 

We choose j = 80, k = tt/2, Jq = 1, a = 6, and 
N = 220. The local spin resides at site 120, and the 
absorbing potential sites are from 180 to 184 with a con- 
stant U = — i at each site. Fig. 3 shows the time evolution 
of both \<fi2,j(t)\ 2 and |</j 3j -(i)| 2 . Part (a) shows the spin- 
up component of the electron. At the initial scattering 
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center (local spin at site 120), a transmitted and small re- 
flected component emerges; the transmitted component 
then propagates toward the absorbing potential where 
some is reflected, very little is transmitted, and most is 
absorbed. Part (b) shows the spin-down component of 
the electron. There is none until the electron interacts 
with the local spin. (Notice the difference in vertical scale 
compared to (a)). Equal amounts are propagated to the 
right and to the left, and eventually most of the packet 
on the right is absorbed by the negative imaginary poten- 
tial. (As in (a), some is reflected and some is transmitted 
as well). 

Each profile in the plot corresponds to a snapshot of (a) 
\ t P2,j{t)\ 2 or (b) \V3j(t)\ 2 , taken at different times from 
t = to t = 65. The wave packet leaves the local spin at 
t = 30 and does not interact with the complex potential 
until t = 40. Eventually the transmitted wave packet 
from the spin-flip scattering interacts with the complex 
potentials; by t — 60 the interaction is complete. No 
interaction occurs beyond t — 60. 

The important dynamics of the local spin can also be 
illustrated by (S z (t)), which indicates the degree of spin 
transfer from the incoming electron to the local spinJ^ 
We plot (S z (t)) as a function of t in Fig. 4 for a variety 
of values of S. In all cases, it is obvious that the first 
increase of (S z (t)} is induced by spin transfer from the 
incoming electron to the local spin as illustrated well in 
Fig. 3b. This is the increase that one will obtain when 
the Hamiltonian has only real- valued potentials. (S z (t)} 
remains unchanged from t — 30 to t = 40, which is also 
consistent with the time evolution of the wave packet. It 
is after t = 40 that an intriguing feature takes place when 
the transmitted wave packet starts interacting with the 
complex potentials. (S z (t)) increases considerably even 
though (i) the scattering is spin-independent, and, more 
importantly, (ii) the wave packet no longer physically 
overlaps with the local spin. This phenomenon could 
be interpreted as a complex potential-driven action-at-a- 
distance. The equation governing (S z (t)) can be derived 
based on Eq. (fTO)) . Since H = H n — iVo, where H 
and Vq is real valued, one can show 



Hi 



i([H ,S z ])-2[(V S z )-(V )(S z )] . (22) 



Since the two interactions do not occur at the same time 
because of their arrangement, we can consider the two 
separately. The first term gives rise to spin transfer 
while the additional increase after t = 40 is attributed 
to the second term. As mentioned earlier, if the ab- 
sorbing potentials were real-valued or zero, there would 
be no additional change in (S z (t)). In Fig. 4, we show 
(S z (t))/S, where S = \ (S Z (0))\, for various values of local 
spin S = 1/2, 1, 3/2, and 2. As shown in the plot, the 
effects of the complex potential are most significant for 
5 = 1/2. On the other hand, for S > 3/2 the effect of the 
potential is less considerable and, when the electron in- 
teracts with the absorbing potential, (S z (t)) of the local 
spin actually decreases. 
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FIG. 3: (Color online) The time evolution of the (a) \tp2,j\ 2 
and (b) \ips.j | 2 components of a wave packet in the presence of 
a spin-flip interaction as well as a complex absorbing poten- 
tial. The initial mean position of the wave packet is xo — 80 
with ko = 7r/2. The lattice size N = 220. In this case, one can 
actually use total angular momentum conservation to argue 
that the other two components are zero. 



Even though the displayed crossover appears for large 
S, a similar crossover can occur even for S = 1/2, as will 
be now explained. The definition of (S z ) (see above Eq. 
(fT!))) ) includes specific combinations of different compo- 
nents of the wave packet. The relevant components here 
are those that strike the absorbing potential region; they 
remain entangled with the local spin, and hence the spin 
is affected too. Let us first examine what happens when 
a wave packet strikes an absorbing potential. Fig. 5 illus- 
trates the flux reduction as a function of impurity poten- 
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FIG. 4: (Color online) The dynamics of (S z ) for S = 
1/2, 1, 3/2, and 2. The local spin is at site 120 (reached at ap- 
proximately time t — 20). The spin- 1/2 electron is coupled to 
the local spin with coupling Jo = 1. The complex potentials 
are at 5 sites from 180 to 184, reached after time 50; their 
strengths are a constant, U = —i. 




FIG. 5: (Color online) Flux reduction for different number of 
complex potential sites. 



tial U, for a varying number of impurity sites (or varying 
width of imaginary impurity potential well). Note that 
we define U = Uq — iVo, and here we use pure imaginary 
potentials only. All imaginary values of the complex po- 
tential are taken to be negative (Vq > 0); positive values 
result in a flux increase, and are of no interest here. 
Analytical expressions can be easily derived for one 



and two sites, by using the expressions in the appendix 
of Ref. [TH for U = — Wq. In the case of a single site, one 
obtains 



Fluxi = \T\ 



1 



Vq sin k 



sinfc 2 + Vbsinfc + (Vb/2) 2 ' 

where k is the wave vector of the plane wave. Note that 
one obtains unit flux both for Vq = and for Vq — > oo; 
otherwise there is a sharp reduction as Vq increases from 
zero, followed by a steady recovery for increasing values. 
This expression is plotted (for k — tt/2) in Fig. 6, and 
is indistinguishable from the result obtained with a wave 
packet with finite width (a — 6). The expression for two- 
impurity sites shows the same characteristics; for k = tt/2 
it is 



Flux, 



4 + 



(2 + 2V + V 2 ) 



2\2 ' 



(24) 



This curve is also indistinguishable from the wave packet 
result in Fig. 6, and results for more absorbing poten- 
tials differ by very little from the case of two impurities. 
Thus it is clear that an optimum imaginary potential 
strength exists to maximize the flux reduction. In Fig. 6, 
in addition, we show the flux reduction for 5 absorbing 
impurities, with and without the local spin, and for two 
different initial conditions for the local spin. In this plot 
the flux is normalized to the flux transmitted through 
the local spin. Finally, in the interest of completeness, 
one can ask whether a real component of the absorbing 
potentials has a strong influence on what has been done 
so far. In Fig. 7 we show the flux reduction as a function 
of the potential U = Uq — iVq] as before only positive val- 
ues of Vq are shown. This result is for a single absorbing 
impurity. Results for more impurities are similar to this 
one. The result for a single impurity is, for k = tt/2, 



Flux = 1 



Wq 



(2 + V ) 



(25) 



and agrees perfectly with the corresponding numerical re- 
sult shown. As Fig. 7 indicates, a large Uq serves to skew 
the Uq = result. For a single impurity the minimum 
flux occurs at Vbmin = \/4 + Uq and the flux reduction 
slowly goes to zero as Uq increases. 

Returning to the problem that includes the spin-flip in- 
teraction, we can define reflection and transmission mag- 
nitudes that leave the local spin in a particular config- 
uration. For example, let Rl and R2 be the reflection 
coefficients with local spin up and down, respectively, 
and let Tl and T2 be the transmission coefficients with 
local spin up and down respectively. With these mag- 
nitudes the expectation value of the z-component of the 
local spin after interaction can be calculated to be 



R1 + T1-R2- T2 



1 



R2 + T2 



2{Rl + Tl + R2 + T2) 2 RI + T1 + R2 + T2' 

(26) 
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FIG. 6: (Color online)Flux reduction at the imaginary po- 
tential, with and without local spin. Expected and numerical 
results for one and two site impurity potentials. 5 sites with 
local spin with different orientations. 




FIG. 8: (Color online) Time evolution of (S z ) of local spin 
with an electron spin up interacting with local spin and 
imaginary potential with different electron-spin interaction 
strength. 




2 4 6 8 10 12 



FIG. 7: (Color online) Flux reduction for complex potential 
with both real and imaginary parts. 



This equation shows that the final value of (S z ) depends 
on a particular combination of these magnitudes. There- 
fore initial conditions and system parameters, which de- 
termine these magnitudes, change the value of (S z ). An 
example of this different behavior arises as a function of 
the Jo coupling of the electron spin to the local spin. The 
expectation value of the z-component of the local spin, 
(S z ), is shown as a function of time in Fig. 8, for a variety 



of values of Jo- As illustrated, the change after t = 50 al- 
ters its characteristic around Jo = 2.3; namely, instead of 
an additional increase, (S z ) decreases when the coupling 
(to the local spin) is sufficiently strong. 

As was the case with varying magnitude of spin (see 
Fig. 4), the local spin can react in a qualitatively very 
different way when the transmitted part of the electron 
reaches the absorbing potential. A summary of the scat- 
tering coefficients after interaction with the local spin is 
displayed in Fig. 9. Note that Jo ~ 2.3 plays a promi- 
nent role; it is the value of coupling for which the spin- 
flip component of the transmitted wave packet peaks as 
a function of Jo- One can readily show that the change 
in (S z (t)) as given by Eq. (|26|) after scattering from an 
imaginary potential (compared with the change before 
scattering) is proportional to T2 — R2. This quantity 
changes sign at Jo « 2.3, and therefore leads to the qual- 
itative change indicated in Fig. 8, independent of the 
value of the absorbing potential. 

The impact of the magnitude of Vq (negative imagi- 
nary absorbing potential) on the magnitude of change is 
illustrated in Fig. \W[ for a couple of local spin initial 
configurations, (a) (S) initially in the -z-direction, and 
(b) (S) initially in the x-direction. 

In figure [TOT a). (S z ) is plotted for a time after the 
transmitted part of the wave packet interacts with the 
imaginary potential, for the case where the initial local 
spin configuration was aligned in the negative z-direction 
(the incoming electron is always spin up) . For low values 
of Jo, the imaginary potential first causes an increase in 
the local z-component spin, followed by a decrease (for 
large enough Vq the local spin is not affected, which is 
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FIG. 9: (Color online) Magnitudes of each component of the 
wave packet with respect to local spin. Rl has the same 
value as Tl, as the created state of electron down, local spin 
up needs to conserve momentum. 



consistent with Fig. 5). As already noted, for larger 
Jo (greater than about Jo w 2.3 in the figure), the ef- 
fect of the absorbing potential is opposite: the local z- 
component spin first decreases and then increases as a 
function of Vq. 

Note that the initial configuration described above is, 
in many ways, quite unique, because the local spin and 
the electron spin are perfectly anti-aligned initially. It 
turns out that for this initial configuration (and for some 
others), J z is conserved, contrary to the general expec- 
tation based on Eq. JTSJ (with A = J z ). In Fig. [TUfb) 
we show an example where J z is not conserved, as one 
would expect in general. To show this more explicitly, in 
Fig. 11 we plot a sequence of results for initial starting 
configurations that sweep through the x-z plane. At the 
two end-points J z remains constant as a function of time. 
For other initial conditions, (J z ) changes as a function of 
time. As Fig. [5] showed, the amount of flux reduction 
due to the absorption is dependent on the value of Vq. 
For large Vq — ► oo, the reduction is diminished by the 
reflection from the potential, and full flux is recovered. 
This effect is visible in (S z ) as well. The main change 
occurs around Vq ~ 1, and full recovery (to the Vq = 
value) occurs for high values of the interaction. 



V. SUMMARY 

Complex absorbing potentials have been utilized as a 
source of absorption in the context of dynamical prob- 
lems (like spin-flip scattering) by investigating the time 
evolution of a particle wave packet under the influence 
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FIG. 10: (Color online) z-component of local spin for different 
interaction strengths of electron-spin and imaginary potential, 
with initial configuration of (a) |+ J.), and (b) | + 0). 



of such potentials. The formulation presented in this pa- 
per is mathematically equivalent to diagonalizing a com- 
plex symmetric matrix because the Hamiltonian we con- 
sider is, in the matrix representation, complex symmetric 
rather than Hermitian. The structure of the formalism 
parallels the one in 'conventional' quantum mechanics 
with Hermitian Hamiltonians. However, the diagonal- 
izability of complex symmetric matrices requires some 
extra attention. Eq. (JTUJ) is central to the results when a 
complex potential is considered. This equation indicates 
that conservation of a physical quantity does not follow 
from the commutability of its operator with the Hamil- 
tonian in the presence of a complex (i.e. non-Hermitian) 
potential. One of the most remarkable examples is the to- 
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FIG. 11: (Color online) (J z ) of the local spin for various local 
spin initial conditions as a function of time. Jo = 1 and 
Vo = 1. 



tal angular momentum J. We showed an example where 
( J z ) is not conserved, and shows significant time depen- 
dence, depending on the initial conditions. 

We also studied the interplay between the spin-flip 
potential and the absorbing potential. The first, while 
Hermitian, nonetheless displays some unusual proper- 
ties connected to entanglement. Essentially, the electron 
wave packet that emerges from the spin flip scattering 
cannot be separated from the local spin state. As a con- 
sequence, when the transmitted portion further interacts 
with an absorbing potential, the state (and expectation 
values) of the local spin is affected. We also studied the 
'absorbing power' of a complex potential as a function 
of its strength (both real and imaginary parts). The re- 
markable 'action-at-a-distance' displayed in Figs. 4 and 
5 is a consequence of an imaginary potential. It would 
be most intriguing if this effect can be observed in a real 
experiment. 

A complex absorbing potential for matter waves can be 
realized by a standing light wave resonant with a tran- 
sition between a metastable state and an excited state, 
which corresponds to an open two-level system. Such a 
scheme has been used in various experiment s 19 ! 20 ' 21 in 
quantum optics. 

Theoretically, complex potentials are used frequently 
in quantum mechanical problems to simulate absorption 
processes. Nonetheless all of their effects including the 
one we illustrate here may not be completely physical. 
However, this is far from clear and requires an experi- 
mental test to delineate the possibilities. As mentioned 
above, the most promising candidate for a physical real- 
ization is a quantum optical system. 
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APPENDIX A: COMPLEX SQUARE- WELL 
POTENTIAL 

For a complex square-well potential, V(x) — Vq for 
|x| < a while V(x) — elsewhere. The general solution 
to the eigenvalue equation, Hip = Etp, can be expressed 
as follows: 

{C cos(kx) + Dsin(kx) if x < —a, 
Fcos(kx) + Gsin(nx) if — a < x < a, (Al) 
Acos(kx) + Bsin(kx) if x > a, 

where A, B, C, D, F and G are complex constants, k = 
^2m(E - Vq) and k = \j2mE. As the Hamilt oman is 
non-Hcrmitian, both k and k lie in the complex plane. 
Since the potential V(x) is even, the Hamiltonian op- 
erator commutes with the parity operator. Hence, all 
non-degenerate eigenstates of the Hamiltonian must also 
be eigenstates of parity. 

Consider first the even states. Evenness of ip implies 
that A = C, D = —B, and G = in Eq. The 
derivative d x ip must be odd, so d x ip\L = —d x ip\-L- How- 
ever, periodic boundary conditions require that d x tp\i, = 
d x tp\-L- Therefore, 

d x ip\ L = k (-Asin(kL) + B cos(kL)) = 0. (A2) 

By the continuity of tp and ipl at x = a, 

F cos(na) + G sin(Ka) = A cos(fca) + B sin(fca) , (A3) 

and 

k (~Fsm(Ka) + Gcos(na)) = k (— Asin(ka) + B cos(fcet)) 

(A4) 

It is easy to search for eigenvalues which satisfy 
cos(fcL) = 0, cos(ko) = 0, or cos(fca) = 0, since the ze- 
ros of these functions are well-known. However, in most 
configurations, the vast majority of eigenvalues do not 
satisfy these conditions, so Eq. (20) must be solved nu- 
merically in the complex plane. This is done by applying 
Newton's method to the following function: 

fn(k) — tan -1 ^— tan(/«a)^ + k(L — a) — nn , (A5) 

for each n. 
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Analysis of the odd states is similar to that of the even 
states. Oddness implies that D = B, C = —A and F = 
in Eq. (|Alj) ■ We apply periodic boundary conditions and 
require continuity of ip and ipl at x = a in order to derive 
the following expression for odd eigenstate momenta k. 
If cos(kL) ^ 0, cos(rea) ^ 0, and cos(fca) ^ 0, we obtain 

fctan(Ka) = — Ktanfc(L — a) . (A6) 

In order to solve this equation in the complex plane, 
Newton's method is applied to the function 

9n{k) = tan -1 ( — tan(«;a) ) + k(L — a) — mr , (A7) 
\k J 

for each n. It can be seen numerically that in most cases, 
f n {k) and g n {k) each have one zero per value of n. Oc- 



where the ± in front implies '+' for even states and '— ' 
for odd states, and jV = {n\n) is a normalization factor. 
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